Iterative Solution of Skew-Symmetric Linear Systems
نویسندگان
چکیده
We offer a systematic study of Krylov subspace methods for solving skew-symmetric linear systems. For the method of conjugate gradients we derive a backward stable block decomposition of skew-symmetric tridiagonal matrices and set search directions that satisfy a special relationship, which we call skew-A-conjugacy. Imposing Galerkin conditions, the resulting scheme is equivalent to the CGNE algorithm, but the derivation does not rely on the normal equations. We also discuss minimum residual algorithms, review recent related work, and show how the iterations are derived. The important question of preconditioning is then addressed. The preconditioned iterations we develop are based on preserving the skew-symmetry, and we introduce an incomplete 2× 2 block LDLT decomposition. A numerical example illustrates the convergence properties of the algorithms and the effectiveness of the preconditioning approach.
منابع مشابه
Convergence Properties of Hermitian and Skew Hermitian Splitting Methods
In this paper we consider the solutions of linear systems of saddle point problems. By using the spectrum of a quadratic matrix polynomial, we study the eigenvalues of the iterative matrix of the Hermitian and skew Hermitian splitting method.
متن کاملA Preconditioner for Generalized Saddle Point Problems
In this paper we consider the solution of linear systems of saddle point type by preconditioned Krylov subspace methods. A preconditioning strategy based on the symmetric/ skew-symmetric splitting of the coefficient matrix is proposed, and some useful properties of the preconditioned matrix are established. The potential of this approach is illustrated by numerical experiments with matrices fro...
متن کاملNumerical Investigation of Krylov Subspace Methods for Solving Non-symmetric Systems of Linear Equations with Dominant Skew-symmetric Part
Numerical investigation of BiCG and GMRES methods for solving non-symmetric linear equation systems with dominant skew-symmetric part has been presented. Numerical experiments were carried out for the linear system arising from a 5-point central difference approximation of the two dimensional convection-diffusion problem with different velocity coefficients and small parameter at the higher der...
متن کاملMinimal Residual Methods for Complex Symmetric, Skew Symmetric, and Skew Hermitian Systems
While there is no lack of efficient Krylov subspace solvers for Hermitian systems, few exist for complex symmetric, skew symmetric, or skew Hermitian systems, which are increasingly important in modern applications including quantum dynamics, electromagnetics, and power systems. For a large, consistent, complex symmetric system, one may apply a non-Hermitian Krylov subspace method disregarding ...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
- SIAM J. Matrix Analysis Applications
دوره 31 شماره
صفحات -
تاریخ انتشار 2009